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O ' Abstract 
(N' 

We investigate the nonlinear instability of a smooth Rayleigh- Taylor steady-state solution (in- 



^. eluding the case of height heavier density with increasing height) to the three-dimensional incom- 

"^ ! pressible nonhomogeneous magnetohydrodynamic equations of zero resistivity in the presence of 

O I a uniform gravitational field. We first analyze the equations obtained from linearization around 
the steady-state solution. Then we construct solutions of the linearized problem that grow in 

^ • time in the Sobolev space H'^, thus leading to a instability result for the linearized problem. 

^ . With the help of the constructed unstable linear solutions and a local well-posedness result of 

^ ! classical solutions to the original nonlinear problem, we can show the instability of the nonlinear 

-)— » 

a 



problem from the dynamics point of view. Moreover, we find that the third component of the 
velocity already induces the nonlinear instability. 

Keywords: Navier-Stokes equations, steady solutions, Rayleigh- Taylor instability, 
incompressible magnetohydrodynamic fluids. 
^ : 2000 MSC: 35Q35, 76D03. 
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^i ' 1. Introduction 

O ' This paper is concerned with the nonlinear Rayleigh- Taylor (RT) instability of a smooth RT 

steady-state solution to the following three-dimensional nonhomogeneous incompressible magne- 
tohydrodynamic (MHD) equations with zero resistivity (i.e. without magnetic diffusivity) under 

on 



the presence of a uniform gravitational field in a domain Q (see, for example, [2|, |2l|, |23|, |25 
the derivation of the equations): 

Pi + V ■ Vp = 0, 

pvi + pv ■ Vv + Vp = (V X M) X M + /iAv - pge^, 

Mt - V X (v X M) = 0, ^ ' ' 

divv = 0, divM = 0. 

Here the unknowns p := p(t,x), v := v(t,x), M := M(t,x) and p := p{t,x) denote the density, 
velocity, magnetic field and pressure of the incompressible fluid, respectively, p > stands 
for the coefficient of shear viscosity, g > for the gravitational constant, e^ = (0, 0, 1) for the 
vertical unit vector, and —ge^ for the gravitational force. In the system (II. ip the equation (ll.ljl i 
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is the continuity equation, while f lLip o describes the balance law of momentum. It is well- know 
that the electromagnetic field is governed by the Maxwell equations. In MHD, the displacement 
current can be neglected [2l|, |23|. As a consequence, f ll.ip -:! is called the induction equation. As 
for the constraint divM = 0, it can be seen just as a restriction on the initial value of M since 
(divM)j = 0. We remark that, the resistivity in (11. ip ^ is zero, which arises in the physics regime 
with negligible electrical resistance, see |3[. In addition, if M = 0, equations fll.ip reduces to a 
standard Navier-Stokes equations of incompressible viscous fluids in the presence of a uniform 
gravitational fleld. 

In this paper we consider the problem of the RT instability in a horizontal periodic domain. 



thus we deflne ^ := (2'7rLT) 



X 



and assume that a smooth RT (steady-state) density profile 



p := pi^X's) E C°°(K) exists which satisfies 

p' G Co°°(M), inf p > 0, (1.2) 

X3GR 

p'{x'l) > for some point X3 G R. (1.3) 

Here ' = d/dx^, and we have written 2itLT for the ID-torus of length 27rL. We refer to |16l . 
Remark 1.1] on the construction of such p. Let M G M^ be a constant magnetic field, then the 
RT density profile p with (v, M)(t, x) = (0, M) defines a steady state to (II. ID . provided (in fact 
the steady pressure p is determined by p) 

V7- - .dp 

Vp = -pge3, I.e., -— = -pg. 

dX3 

We point out that by virtue of the condition (II. 3p . there is at least a region in which the RT 
density profile has larger density with increasing X3 (height), thus this will lead to the classical 
RT instability as shown in Theorem 11.11 below. 
Now, we denote the perturbation by 

- 0, N = M - M, 



u 



q = p-p, 



/iAu+(V X N) X (N + M), 



;i-4) 



g = p- P 
then, {g, u, q) satisfies the perturbed equations 

ft + u ■ V(^ + p) = 0, 
{g + p)ut + {g + p)u ■Vu + Vq + gge^ 
Ni = Vx (ux (N + M)), 
^ divu = 0, div N = 0. 

To complete the statement of the perturbed problem, we specify the initial and boundary condi- 
tions: 

(ftu,N)|t=o = (^o,uo,No) infi (1.5) 

and 

lim (u,N)(t,x',a;3) = for any t > 0, (1.6) 

|X3|— >-+00 

where we have written x' = XiGi + ^262, ei := (1,0,0) and 62 := (0, 1,0). Moreover, the initial 
data should satisfy the compatibility conditions divuo = and div No = 0. 

If we linearize the equations (ll.4p around the steady state (p, 0), then the resulting linearized 
equations reads as 

gt + p'us = 0, 

put + Vg + 5-^63 = pAu + (V X N) X M, 

Nt = V X (ux M), 

divu = 0, divN = 0. 



;i-7) 



The RT instability is well known as gravity-driven instability in fluids when heavy fluid is 
on top of light one. The linear instability for an incompressible fluid was flrst introduced by 



Rayleigh in 1883 [28[. The analogue of the RT instability arises when the fluids are electrically 
conducting and a magnetic fleld is present, and the growth of the instability will be influenced 
by the magnetic fleld due to the generated electromagnetic induction and the Lorentz force. 
Some authors have extended the partial results concerning the RT instability to the case of 
MHD fluids by overcoming difficulties induced by presence of the magnetic fleld. For instance, 
Kruskal and Schwarzchild in 1954 flrst showed that a horizontal magnetic fleld has no effect on 
the development of the linear RT instability ^]. Then the influence of a vertical magnetic fleld 



was investigated by Hide in [13[ where the effect of flnite viscosity and resistivity was included 
and his analysis was encumbered with many parameters. By a variational approach, Hwang in 
2008 studied the nonlinear RT instability of f ll.4p -f lT3]) for the inviscid case (i.e. /i = 0) in a 
two-dimensional periodic domain [14!]. 

To our best knowledge, however, it is still open mathematically whether there exists a unstable 
solution to the nonlinear RT problem fll.4p -f fTT6|) of 3D viscous MHD fluids. The aim of this article 
is to prove rigorously the instability for the nonlinear RT problem f ll.4p - fll.6p . The main result 
reads as follows. 

Theorem 1.1. Let the critical number 

Mr:= 



\ 



sup f. > 0. (1.^ 

^6-ffi(K) J]j \W I 0.x 



Assume that the RT density profile p satisfies /il.^) and /il.3\) . and 

{ Mei, M 7^ 0, constant 

t Meg, |M| e (0,Mc), constant. 

Then, the steady state (p, 0, M) of {l.^^ - ^TT^ is unstable, that is, there exist positive constants 
A*, £0 o,nd L, and a triple (^o,Uo,No) G H°°{Q,) := n^Qif'^(f2), such that for any 5 G (0, t) 
and the initial data {qq, uq, Nq) := {Sgo, 5uo, SNq) there is a unique classical solution {g, u, N) of 
( [i.y^[ j- pT^) on [0,T™^^), but the third component of the velocity satisfies 

\\u3(T^)\\L^(n) > eo for some escape time T^ := — In^^ G (0,r"'''^), (1.10) 

A* 

where T™^^ denotes the maximal time of existence of the solution (f), u, N). 

Remark 1.1. We should point out that since the above Sobolev spaces H^{Vt) are deflned on 
the horizontal periodic domain Vt = (27rLT)^ x M, the solution {g, u, N) constructed in Theorem 
11.11 is horizontally periodic. One should keep this in mind in what follows. 

Throughout the rest of this article we shall repeatedly use the abbreviations: 

lYm,p ._ lYm,p^^^^ jjm ._ ffm^^^^ jjoo ._ n^^^iy™(fi), L^ := L^ (Q) , 

II • ||vK™.p •= II ■ ||iy™.p(n)! II ■ ||h™ •= II ■ ||_f/™(n)i II ■ IIlp '■= II ■ ||lp(q)5 etc. 

Remark 1.2. Theorem 11.11 gives the instability of the third component of the velocity. We 
mention that the instability of the third component of magnetic fleld (i.e., N^) can be obtained 
for the case of M = Me^ thanks to Remark |2.4[ More precisely, one can use the linear instability 
of A^3 in Remark 12.41 to show the nonlinear instability of A^3 in the original problem in the same 
way as one utilizes (I2.60p to get the nonlinear instability (ll.lOp . 



Remark 1.3. Theorem 11.11 also holds for the general horizontal magnetic field M = (Mi, M2, 0). 
In fact, rotating the o-xy coordinates properly so that M = (M, 0, 0) with M = \/ Ml + Mf , 
we have the same case as in Theorem 11.11 under the rotated coordinates, since the L^-norm of 
the third component of velocity is invariant under the horizontal rotation. In addition. Theorem 
11.11 also holds when M = N = (the case of incompressible viscous fluids), hence the weak RT 



instability for incompressible viscous fluids given in [16| can be further shown to be the strong 



RT instability as (ll.lOp in the case of the horizontal periodic domain. 

Remark 1.4. The number Mc in ( 11. Sp may be inflnite. In fact, we can construct a RT density 
proflle p satisfying (11.21) . (11.31) and Jjgp'dx > 0. Let supp' C (—1,1) with / > 0, and 



Then V"™ (a;) G H\R) and 



1, X e [—/,/], 

1 + (x + l)/n, X G (— / — n, — /), 

1 — (x — l)/n, X E {1,1 + n), 

X ^ (— / — n,l + n). 






> gn p'dx — )■ 00 as 77. 

Jr 



—7- CXD, 



which implies M^ = 00. Thus we see that the vertical magnetic fleld in this case has no influ- 
ence on the development of the nonlinear RT instability, provided p(+oo) := \im.x.^^+oo pi^s) > 
p(— 00) := lim2^3^_oop(x3). Moreover, in such a case, we can infer that the critical frequency 
constant \^\^^ deflned by (12.151) is equal to zero. 

Remark 1.5. In [7| Hwang and Guo proved the nonlinear RT instability for two-dimensional 
nonhomogeneous incompressible inviscid flows (i.e. p = and M = N = in (II. 4p . (II. 5p ) with 
boundary condition u ■ n|aQ/ = where Q' = {(xi,X2) G M^ | — I < X2 < m} and n denotes the 
outer normal vector to dQ'. Our instability result Theorem 11.11 differs from that of Hwang and 
Guo in that only HMsHia is needed here to describe the instability. This is also different from that 
of Hwang [7| , in which the instability for an inhomogeneous incompressible inviscid MHD fluid 
is described by the norm ||(f?, u, N)||j;^2(q/). 

The proof of Theorem 11.11 is divided into four steps given in Sections [2H5l (i) First, we 
make an ansatz to seek the "normal mode" solutions of the linearized equations (II. 7p . which are 
supposed to be grow exponentially in time by the factor e^'-^-'* with ^ G M^ being the horizontal 
spatial frequency and A(.^) > 0. This reduces the equations to a system of ordinary differential 
equations (ODEs) deflned on M with X{^) > for some ^ (see (12. 4p . (12. 5p ). Those points ^ 
constitute a solvable domain. Because of presence of the magnetic fleld, the solvable domain is 
not a ball for the horizontal case as in [10|, resulting in additional technical difficulties. In order 
to circumvent such difficulties, we shall introduce new notion of the critical frequency function 
S{^) and critical frequency constant |^|*^ (see (I2.14p . (I2.15P ) to define the solvable domain A^. 
Thus, by careful constructing the solutions and adapting the modified variational method in [lO| , 
we can also solve the ODEs for any given ^ G A^ and obtain a normal mode with A(^) > 0, 
providing thus a mechanism for the global linear RT instability. Using the normal modes, we can 
further construct real- valued solutions of the linearized Cauchy problem (ll.5p - (ll.7p that grow in 
time, when measured in H'^[Q) for any k > 0. In particular, the third component of velocity 
in the linear solutions are non-zero functions, which plays a key role in the nonlinear instability 



in f ll.lUp . (ii) In Section [31 we further show that A is indeed a sharp exponential growth rate 
for general solutions to the linearized problem. The main idea in the proof was first used by 



Hwang and Guo for the inviscid incompressible fluids, and then adapted by Guo and Tice [10 
for viscous fluids. Thus, based on Duhamel's principle, we can bound the nonlinear terms by the 
sharp exponential growth rates (see (15. 7p ). (iii) In Section HI we derive the integrand form of 
Gronwall's inequality of high-order energy estimate S{g,u, N) for the perturbed problem (11.41) - 
(11.61) . which makes the escape time occuring before the blow-up time of the classical solutions. 
Since the equilibrium state of the magnetic field M is a no-zero vector, we shall introduce a 
new technique to deal with the terms including M in the energy estimates, see Subsection 14.41 
(iv) Finally, in Section |5l with the help of the results established in Sections [JHU we modify the 
approach in [^ to establish the instability of the nonlinear problem in the sense of (ll.lOp . We 



mention that the approach in [8[ has been also used to treat the instability of other problems. 



see [g, [151, |31[ for examples. 

We end this section by briefly reviewing some of the previous results on the nonlinear RT 
instability for two layer incompressible fluids separated by a free interface (stratified fluids), 
where the RT steady-state solution is a denser fluid lying above a lighter one separated by a free 
interface. When the densities of two layer fluids are two constants, Wang and Tice [31[ proved 
the (local) existence of nonlinear unstable solutions in a horizontal periodic domain T^ x (—6, 1), 
where the instability term is described by the sum of L^-norm of the velocity and the moving 
internal interface. Priiss and Simonett used the C°-semigroup theory and the Henry instability 
theorem to show the existence of nonlinear unstable solutions in the Sobolev-Slobodeckii spaces 
in I 



27| . where the instability term is described by the sum of ||u||^^2-2/p.jj3. and ||h||^3-2/p 



(see [271, Theorem 1.2] for details). When densities of two layer fluids are variable, to our best 
knowledge, the (local) existence of solutions to the nonlinear problem (ll.4p - (ll.6p still remains 
open, consequently, the strong nonlinear instability is still open. For compressible fluids there are 
very few results on the nonlinear RT instability. Guo and Tice proved the instability of immiscible 
compressible inviscid fluids in the frame of Lagrangian coordinates under the assumption of the 



existence of solutions [12| . This is in some sense a compressible analogue to the local ill-posedness 



ipre 

i. 



of the RT problem for incompressible fluids given in 

Finally we mention some results on the instability for stratified MHD fluids. Wang [30 
introduced the critical number of stratified MHD fluids (denoted by M^) to investigate the linear 
RT instability of stratified MHD fluids in a infinite slab domain M^ x [—1,1). Later, Jiang et 



al. [17[ further showed the weak nonlinear RT instability of stratified MHD fluids and found 
that M^ = ^ygl(g_^_ — Q-)/2, where ^+ > Q-, and ^+ resp. ^-_ denotes the density of the upper- 
resp. lower-layer fluid. Obviously, M^ — > oo as the height 2/ — )■ oo. Hence the critical number 
of stratified MHD fluids in M^ is infinite. This means that the vertical magnetic field may have 
no influence on the development of the RT instability in M?. Due to some technical difficulties, 
however, we do not know whether Theorem 11.11 still holds for the corresponding stratified MHD 
fiuids. 

2. Construction of solutions to the linearized problem 

We wish to construct a solution to the linearized equations (11.70 that has growing if '^-norm 



for any k. We will construct such solutions via some synthesis as in [10[ by first constructing a 
growing mode for any but fixed spatial frequency. Moreover, we shall introduce the techniques 
of the critical frequency function and critical frequency constant to carefully analyze the terms 
involving with the magnetic field. 



2.1. Linear growing modes 

To start with, we make a growing mode ansatz of solutions, i.e., 

g{t,x.) = p(x)e^*, u(t,x) = v(x)e^*, g(t,x) = p{x.)e^\ N(t,x) = M(x)e^* for some A > 0. 

Substituting this ansatz into (II. 7p . and then ehminating p and M by using the first and third 
equations, we arrive at the time-independent system for v = {vi,V2,V3) and p: 



with 



A^pv + AVp = (V X (V X (v X M))) x M + A/iAv + gpv-^e^, 
div V = 



hm (v,M)(x',a;3) = 0, 

|a:3|— >-+co 



(2.1) 



(2.2) 



where M is given by (II. 9p . After a straightforward calculation, we find that 
V X (V X (v X M)) X M = V X (M ■ Vv) x M 



M^ 



(0, dlu2 - dl^ui, df^u-s - dl^ui), for M = Me 



1) 



(2.3) 



{d'^ui - di^us, dl^U2 - dl^uz, 0), for M = Me-^. 
We fix a spatial frequency ^ = (^1,^2) G I^^, and define the new unknowns 

{}i(x) = -i^(x3)e''''-«, 52(x) = -i^(x3)e^^' -^ U^) = ^^{x.y'^'^ , p{^) = Ttix.y'^'^. 
Then, in view of (I2.ip - (I2.3I) . we see that ip, 6, ip and A satisfy the following system of ODEs: 



f A^py. - A^ivr + Ap(|e| V - V") = M'^Bu 
X^pe - A^vr + Ap(|eP^ - e") = M^Sa, 
A^pV^ + Att' + ApdePV^ - ^") - gp'iJ = M^Bs, 



(2.4) 



with 



where 



ip{—oo) = 9{—oo) = ip{—oo) = Lp{+oo) 



+00 



ip{+oo) 



0, 



(2.5) 



( B,\ 

B2 

\bJ 



\ 



for M = Mei, 



for M = Me 



3- 



/ 

I (/' + ilv) \ 

r + 6^') 


Eliminating vr from the third equation in (12. 4p . we obtain the following ODE for ip 

^' ^' ^'^'^^ + ''l(^""-|elV'), forM^Mes, 



V 



/ 



=Ap(ieiv-2ieiV' + v^" 



(2.6) 



6 



with 



ip{-oo) = 7/''(-oo) = ip{+oo) = ij'i+oo) = 0. (2.7) 

Next we use the modified variational method to construct a solution of fl2.6p . (12. 7|) . This 



idea can be found in Guo and Tice's paper on compressible viscous stratified fiows llOj, and is 
adapted by other authors to study the instability for other fiuid models j^, |l5|, |l8|, |30|. To this 
end, we now fix a non-zero vector ^ G M^ and s > 0. From (I2.6p and (12 .Zp we get a family of the 
modified problems 






{21 



(t/;""- I^IV'), forM = Me3, 
coupled with the condition fl2.7l) . We define the energy functional of (12.81) by 

E(^) = |epi?o(^) + 5i?i(^) (2.9) 

with an associated admissible set 



where 



A= \i)eH'^{M) 



Jr^i 



j{i^)--= / p(ieri^r + iv^r)dx3 



(2.10) 



^o(V') = M^ 



'''"'f + ^ 



km+ 



\r 



i?p 



- fifpV dx3, for M = Mei, 
^fpVMxs, for M = Mes, 



f2.111 



(2.12) 



Ei(v^) = /i / (4ien^r + iiei V + ^'T)dx3. 

^R 

Thus we can find a — A^ (depending on ^) by minimizing 

-A2(0 = «(0:=mf ^(V^). (2.13) 

■ip&A 

In order to emphasize the dependence on s G (0, oo) we will sometimes write 

ENj.s) := ENj) and ais) := inf ENj.s) < +cx). 

Before constructing the growth solutions, we shall introduce some preliminary results, which 
will be used in Subsections 12. 2f[2^ Let the critical frequency function S{^) for the horizontal case 
and the critical frequency constant |^|^^ for the vertical case be given by the following variational 
forms 



^(0 : = 
forO< \Mii\/\i\ < Me, and 



sup 



^|eP^p'^Ma;3/(M^0^-/Ml^fdx3 



Jr 



pdxa 



(2.14) 



I S I vc 



7 



for \M\ G (0,Me 



(2.15) 



respectively, where M^ is given by fll.Sp . and 

Now we define a solvable domain for a growing solution 

r {eeM^ I |eiM|/|e|G(0,M,), |el<5(0}u{ei = 0}\{0} forM = Mei, 
A*^ = < _ (2.16) 

I {^ e M^ I \i\^i < lel} for M = Meg with \M\ e (0,M,). 

By virtue of the definition of M^, it is easy to verify that the above two definitions fl2.14p . (12.151) 
make sense. Moreover, |^|^ and S{^) are finite. Next we introduce some properties concerning 
with the critical frequency and the solvable domain A^. 



Proposition 2.1. The supremum in ^2.14 ) is achieved for each ^ and M with S,iM ^ 0. More- 



over, S{^) is continuous on D := {^ G M | < |M^i|/|,^| < M,,} for any given M ^ 0, and 

S{^) ^ +00 as j^^O and ^eB. (2.17) 

Proof. We rewrite (12.14p as 

5^(0 = sup g(7A), 

-4^,2 



where 



QW: ''^' 



(li^//'^'^^--/j*'|^^- -^-^{^'' 



i2(K) = 1}- (2.18) 



Then, it is easy to see that S'^{^) > 0, since < |M^i|/|^| < Mc. Let ipn G Al^ be a minimizing 
sequence of S'^{^), i.e., lirasup^^^ Q{ipn) = sup^ ^ QW^ "we have from (I2.18P that ■^/'n is uni- 
formly bounded in H^{M.) on n. Hence there exists a ipQ & H^{M.), such that tpn -^ V'o weakly in 
if^(M) and strongly in LfQ^(M). Hence, 

< supQ{ij) = limsupQ(V^n) < Q(^o) and < ||V^o||l2(r) < liminf ||^„||i2(iR,) = 1. 

Aj^2 n— >cjo n—KX) 

We proceed to verify that ||'?Ao||l2(ir) = 1. Suppose by contradiction that ||'?Ao||l2(ir) < 1, then we 
may scale up V'o by a > 1 so that aipo G Al'^. From this we deduce that 

supQiip) > Qiaipo) = a'^Qiipo) > a^supQiip) > supQiip), 

which is a contradiction. Hence ||V^||l2(m) = 1, which shows that Q{il)) achieves its infinimum on 
Al^. 

Next, we prove the continuity of 5'^(^o) for each given ^o in '^- Letting ,^ G D — ?► ,^0) we have 
I^P/'Ci ~^ I'CoP/^oi' where ^oi represents the first component of ,^o- Without loss of generality, 
we assume that |CP/Ci and |Co|^/Coi belong to a finite interval (a, 6) with a > 0. Denote 5 := 

leive? - leoiveo'i, then ^ ^ o as e ^ ^o. 

On the other hand, Q{ip) achieves its infinimum on Al^, i.e. for any ^ G D, there is ip^ G Al^, 
such that 

91^1 f 7,Li.2a^ f l„;.'|2. 



S\0 ■■= (^ /_ pWx, - I l^^l^dxa. (2.19) 



Substituting |^o| V-^oi = \^?/^i - ^ into ([211]), one has 
Similar to fl2.20p . we obtain 



(2.20) 



S\^o) 



9\^\^ f -J.,.2 ^. f L/./ ,2^. 9^ 



{Mi, 






which, together with (I2.20p . yields that 
Hence, 



-2(t^ c2^^M / ^^ II ./I 



l^nO - S\i,)\ < ^IIp'IIlo^(m) ^ as <5 ^ 0, 

and S'^{i) is continuous at each given point ,^o ^ ID- Finally, (I2.17P obviously holds by the 
definition (12 .14^ . The completes the proof. D 

Proposition 2.2. Let M he defined by 112. 16\) . then 

(1) the set A^ is symmetric on x-axis and y-axis in M^, respectively; 

(2) the set M is a nonempty open set in M^; 

(3) AS n [L-^Zf ^ 0. 

Proof. It is easy to see that the above assertions obviously hold. We remark that we have used 
(I2.17P and the continuity of S'(^) in the verification of the second assertion for the horizontal 
case. n 

Proposition 2.3. Let H^ 0. If ^ E A^, then there exists a 

H^(M) for the horizontal case, 



^ A^vc for the vertical case, 

such that either Eq^iJjq) < 0, or l^pE'olV^) ^ for any 

H^(R) for the horizontal case; 

t2i 



1^ i/^(M) for the vertical case. 

Proof. The above assertions in fact follow from the definitions (II. Sp . (12. lip and (I2.14p -( |2TT6|1 . 
here we only give the proof of the horizontal case for the reader's convenience. 
Let ^ G A^. If ^1 = 0, then obviously, there exists a ipQ, such that 

E, = j^ [(M6)'(|^oP + ^) -^pVo']dx3 = j^-g-p'^ldx, < 0. 



If ^ satisfies |M^i|/|^| e (0, MJ and |^| < ^(0, then by virtue of the definition of (ETTD . there 
also exists a ipo, such that 



En 



(M6)Ml^oP+''^°''^ --'^''^ 



lei^ 



-gp% 



dxs < 0. 



Summing up the above discussions, we see that there exists a ■j/'o, such that i?o(V^o) < for ^ G A^. 
Let i ^ Asu{0}, then ^ can be divided by two cases: (i) |M{i|/|{| e (0,Mc) with |^| > S(0, 
and (ii) |M^i|/|^| > M^ if M^ < oo (noting that this case will not appear if M^ = oo). For the 
first case, in view of the definition of f l2.14p . we find that 



\^\'EoW = \^N |(Mei)ni^P + ^ 



gp'r 



dx3 > for any ip G H^( 



Finally, for the second case, by the definition of fll.Sp . we have 

[ {^^^\f\'-9p'^')dxs>0 for any ^Gifi 
which yields 



ler^o = lep 



(M^i 



p+'^ 



/|2 



lei^ 



- ^pV^ 



dx3 > for any ip G H^( 



Summarizing the above discussions, we see that \^\'^Eo{iIj) > for any ip G H^(M.) ii ^ ^ A^UJO}. 
This completes the proof for the horizontal case. D 

2.2. Solutions to the variational problem 

In this seusection we show that a minimizer of (12.131) exists for the case of inf_4 £'(■?/', s) < 
which will be shown to be true for sufficiently small s in Proposition 12.51 below, and that the 
corresponding Euler-Lagrange equations are equivalent to (12.71) . (12.81) . 

Proposition 2.4. For any fixed s > and C, with \^\ ^ 0, we have 

(1) inf^e^E(V^, s) > -oo. 

(2) if there exists a ip E A, such that E{iIj) < 0, then -E(^) achieves its infinimum on A. 

(3) let ip be a minimizer and — A^ := E{ip), then the pair (ip, X^) satisfies \2.1\) , h2.8\) . More- 
over, ^ G i7°°(M) := n^oi7'=(M). 

Proof. We only show the proposition for the horizontal case, the vertical case can be dealt with 
in the same manner. 

(1) Noticing that for any ip E A, 



Ei^) > -g\^\' [ p'^p'dxs >-9 ^ [ Piei V'dx3 > -g 

JR P L^(M.) JVL 



(2.21) 



we see that E is bounded from below on A by virtue of (II. 2p . This proves (1). 

(2) We proceed to show (2). Let ipn E A he a. minimizing sequence, then E{%pn) is bounded. 
This together with (12. 9p and (I2.2ip implies that %pn is bounded in H'^{W). So, there exists a 
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ipo G H^(R), such that ipn -^ i'o weakly in H'^{R) and strongly in Hl^^iM.). Moreover, by the 
lower semi-continuity, fll.2p and the assumption that E{ip) < for some V^ G ^, we deduce that 

E{ipo) < \immiE{ipn) = inf E{iIj) < 0, and < J(V^o) < 1- 

Suppose by contradiction that J{ipo) < 1. By the homogeneity of J we may find an a > 1 so 
that J^aipo) = 1, i.e., we may scale up tpo so that aipo G A. From this we deduce that 

E{aiJo) = a'^Eitpo) < a^inf E < inf E < 0, 

which is a contradiction since aipQ G A. Hence Jiipo) = 1 and ipQ G A. This shows that E{iIj) 
achieves its infinimum on A. 

(3) Finally we prove (3). Notice that since E{iIj) and Jiip) are homogeneous of degree 2, 
fl2.13p is equivalent to 

a{s)= inf ^. (2.22) 

For any r G M and ^/J G H'^{M.) we take ^(r) := -ijj + rip , then (12.221) implies 

E(^(r)) + AV(^(r)) >0. 

If we set I{t) = E{iIj{t)) + AV(^(r)), then /(r) > for all r G M and /(O) = 0. This implies 
/'(O) = 0. By virtue of (12.91) and (12.101) . a direct computation leads to 



-^K (2.23) 

=^ier [ pi^Hx3-x' /p(ierv^^+^v)dx3. 



where we have used the upper boundedness of p. 

By further assuming that ip is compactly supported in M, we find that ip satisfies the equation 
(12.81) in the weak sense on R for the horizontal case. In order to improve the regularity of ■?/', we 
rewrite (12.231) as 

rrdxs =- [ (gi^i'p'i, - x'm'pi' - {p^y) 

sp Jr V 

+spim'r - leiV) + (M^inr - lerv^)) hxs (2.24) 



fipdxs. 

For any n > 1, let V'l.n, i^2 G C^(M) satisfy ipi^nixs) = 1 for |x3| < n. If one takes 
ip = 4'i,n /T^ ^2dy in (I2.24p . then one has 

{^,J")^',dx, = / /V^i,„ / ^2d2/ - ^''r / ^2di/ - 2^[J"^2 dX3 



^ — oo 



iflpl,n - ij'lni'")dy - '2'lpln'^" ] V^2dX3, 
X3 
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which, recalhng ip G H'^{M.), imphes ip" G iJiQ^(]R) and 



+00 



ip'" = {i^i,nip"y = / (/V'i,n - i^in'^")^y fo^ ^ny X3 with Ixsl < n. 



X3 



Integrating by parts, we can rewrite fl2.24p as 

<v'dx3 =- [ (^leppv - x'm'p^ - ip^' 

sp Jk V 



+./i(2|e| V - lelV) + (M6)'« - lel V) ) ^da:3 



which, keeping in mind that ^j G H'^(R), yields V^"" G L'^(R). Hence V^ G i7f„,(M) nCf^'y^(M), and 
ip (00) = -i/^ (00) = tp (00) = 0. Using these facts. Holder's inequality, and integration by parts, 
we conclude that 

llV^'lli^W = / l^'"l'dX3 = - [rr"dx, < ||^"||L^(M)liriU^(R), 

Jr Jr 

i.e., ijj'" G L2(]R). Consequently, V^ G if^(R) solves (ET]), (IMD- This immediately gives ip G 



Next, we want to show that there is a fixed point such that X = s. To this end, we first give 
some properties of a{s) as a function of s > 0. 

Proposition 2.5. The function a{s) defined on (0, oo) enjoys the following properties: 

(1) a{s) G C[q^(0, oo) is nondecreasing . 

(2) for any C, G A^, there exist constants ci, C2 > depending on g, M , p, fj,, and C,, such that 

a{s)<-ci + sc2. (2.25) 

Proof. We still give the proof for the horizontal case only, and the vertical case can be dealt 
with in the same way. 

(1) Let {V'^j} C ^ be a minimizing sequence of inf^g_4 £'(■?/', S2). From (12. 9p and (I2.12p it 
follows that 

«(si) < limsupi?('?/'^2 5 "^i) — linisup-E'('?/^"2) "^2) = «(s2) for any < Si < S2 < 00. 

n— >oo n— >oo 

Hence a{s) is nondecreasing on (0, 00). Next we shall use this fact to show the continuity of a{s). 
Let / := [a,b] C M"^ be a bounded interval. In view of (I2.2ip and the monotonicity of a{s), 
we know that 

\a{s)\ < max < !«(&)!, (? ||pVpIIl°=(R) f < 00 for any s G I. (2.26) 

On the other hand, for any s G /, there exists a minimizing sequence {ip'g} C ^ of inf^g_4 E{tp,s), 
such that 

\a{s)-Eii,:,s)\<l. 

Making use of fl2:9l)- fl212l) and JKm . we infer that 

S S ./lij S ./id 



^ l + max{\a{b)\,g\\p'/p\\^^^^^} ^ 
~ a a 
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:= K. 



For Sj G / (i = 1, 2), we further find that 



a{si) <\imsnp E{(p'^^,si) <\imsnp E{iIj^^, S2) + \si - S2\limsnp Ei{iIj^J 

<a{s2) + K\si - S2I. 

Reversing the role of the indices 1 and 2 in the derivation of the inequahty f l2.27p . we obtain the 
same boundedness with the indices switched. Therefore, we have 

\a{si) - a{s2)\ < K\si - S2I, 

which yields a{s) E C^^^{0, 00). 

(2) Since ^ G A^, by virtue of Proposition 12. 3[ there exists a ip E H^(R), such that 

Eoii^) = jj^M^iY (l^r + ^) - gp'^'dxs < 0. (2.28) 

On the other hand, H'^{M.) is dense in H^(R), thus there is a function sequence ■?/'„ G H'^{M.), so 
that 

^jJn^^ stronly in if^(R). (2.29) 

Putting (12.281) and (12.291) together, we see that there is a subsequence ipno G H'^iM.), such 
that 'ipno ^ and 

^o(^no) = / [(Mei)^(|^.or + %^) - 9M,]dx. < 0. 



Thus, we have 



E(^) ^ E(^„J ,„2^o(^no) , ^l(V^no. 



«(s) = inf E{ij) = inf ^^ < ^^^P^ = ler 7-"°/ + 5 V- := -Ci + SC2 

for two positive constants Ci := ci{g,M,p, |^|) and C2 := C2{g,M, fi,p, |^|). This completes the 
proof for the vertical case. D 

Given ^ G A^, by virtue of (12. 25 p . there exists a sq > depending on the quantities g, M, p, 
p and ^, such that for any s G (0, Sq], a(s) < 0. Let 

6^ := sup{s I a(r) < for any r G (0, s)}, (2.30) 

then (3g > 0. This allows us to define A(s) = ^J—a{s) > for any s E S^ := (0, 6^). Therefore, 
as a result of Proposition 12.41 we have the following existence for the modified problem (12. 7p . 
(O. 

Proposition 2.6. For each $, E A^ and s E S^ there is a solution i/j = iIj{^, x^) ^ with 
X = A(^, s) > to the problem ^2^, HJj. Moreover, ijj E H°°{R) n A. 

Now, we can use Proposition 12. 5^ (I2.2ip and (I2.30p to find that A(s) G C^^^{S^) is nonincreas- 
ing, A(s) < ^/g\\^/p'/p\\L°°{R), lim5^oA(s) > and lim^^e^ A(s) = if &\^\ < +00. Hence, we 
can employ a fixed-point argument to find s G 5^ so that s = X{^,s), and thus obtain a solution 
to the original problem (12. 6p . (12. 7p . 
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Proposition 2.7. Let ^ E A^ , then there exists a unique s E S^, such that A(^, s) = ^^/~a(s) > 
and s = A(^, s). 



Proof. We refer to 10|, Theorem 3.8] (or 30|, Lemma 3.7]) for a proof. 



D 



Moreover, in view of Propositions 12 . 6 1 and 12 . 7[ we conclude the following existence for the problem 
(El-dST]). 

Theorem 2.1. For each ^ E A^ , there exist ip = "^(^5 2:3) ^0 and A(^) > satisfying Ii2.6\) and 
( f^ . Moreover, ip E H°°(R) n A. 



We end this subsection by giving additional properties of the solutions established in Theorem 
12.11 in terms of A(^), which show that A is a bounded, continuous function of ^. 



Proposition 2.8. The positive function A : A^ — t- M+ is continuous and satisfies 



sup m < J 9 \\p'/p\\ 



(2.31; 



Proof. The boundedness of A in fl2.3ip follows from fl2.2ip . As for the proof of the continuity 
of A, we still give the proof for the horizontal case only. 

First, let ^0 ^ -^ be arbitrary but fixed, and ^ — )• ^o- Without loss of generality, assume 
^ E A^, since A^ is an open set by Proposition 12.21 Then there exists an interval [a, b] C M^ so 
that 1^1 and |^o| e (a, 6). Let 6 = |^p - |^oP, then 5 -^ as |^| ^ |^o|- 

(i) We begin with the proof of the following conclusion: 

lim a{^, s) = a{C,o, s) for any s E S^. (2.32) 

By virtue of Proposition 12. 6[ for any C, E A^, there exists a function ip^ E A, such that 



«(0 



^M4ien^[P+iieP^c+^[f)+(Mei)^^"^'2'-'2 



+ I^^I')-^I^I'pV| dx3<0, (2.33) 



which, together with (I2.10p . yields 



\mm < C3, 



where C3 depends on g, M, /i, p, a, b and s. 

To deal with the term involved with |^i|, we denote 61 
|^|2 = |^o|2 + § and l^ip = l^oiP + Si into fl03D resuhs in 



(2.34) 



16 P ~ I'CoiP- Substitution of 



«(0 = / sp(4ieoi V^P + m% + ^i\') + {Muf (leoi Vd' + \^\?) - 9\io?p^. 



-1 1 2 



+ 5f{5, i,^) + (5i/i(^5) > a{i,) + 5f{5, ^5) + ^1/1(^5), 



dX3 



(2.35) 



where 



/(5,^«) 



and 



s/x(4|^' |2 + 25^g(|eoP^5 + i^'l) + 5^1 + (Meoi) 



M5,^5)= M'{{\U' + 5)\^P^\' + W^\')<lx,. 



' - g-p'^l 



dX3 
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By Holder's inequality and f l2.34p . we can bound 

\f{6,ilj^)\ + \h{6,ip^)\ < C4 for some constant C4. (2.36) 

Similar to f l2.35p and f l2.36p . we also have 

«(eo) > «(0 - 5fi-5,ij^,) - 5M^^,) (2.37) 

and 

\f(-S,ijJ\ + \h{-S,^^,)\<c,. (2.38) 

Combining (105]) with (jOZD, we get 

6fi-S,i;J + S^h{i;J > a(0 - ^(^0) > Sf{5,i;^) + 5MA), 
which, together with (12.361) and (12.381) . implies (I2.32p . Hence, 

lim A(^, s) = X{C,o, s) for any s G S^, (2.39) 

because of A(^, s) = a/— a(^, s). 

(ii) In view of (12.390 and Proposition 12.71 we see that for any £ > 0, there exists a 5 > 
such that |A(^,S5J - A(^o,S5o)| < ^ and s^^ = A(^o,S5o) = V-«(^o,s?o) ^^r any ^ G 5|^ := 
{y G M^ I |y — ^o| < S} C. A^. On the other hand, for each |^| > 0, A(s) is nonincreasing and 
continuous on S^, and there exists a unique s^ G S^ satisfying A(,^, s^) = sg > by Proposition l2.7[ 
Consequently, we immediately infer that |A(^, s^) — X{^o, sgo)| < e with s^ = A(^, s^). Therefore, 
A(^) is continuous. This completes the proof of the proposition. D 

2.3. Construction of a solution to the system ^2.4^ , Ii2.5\) 

A solution to (12. 6p . (12.70 gives rise to a solution of the ODEs (12.40 . (12. 5 p for the growing 
mode velocity u as well. 

Theorem 2.2. For each ^ G A^, there exists a solution {(f,6,ip,TT) := (</), 6',-?/', 7r)(^,X3) with 
A = A(0 > to [K4\ ), [K^\) . and the solution belongs to H°°{M). Moreover i}) ^ A. 

Proof. We still give the proof for the horizonal case only. First, in view of Theorem 12. 11 we have 
a solution {ip, A) := (?/'(^, Xg), A(0) satisfying ([2SD, (EZD- Moreover, A > and ?/> G Ar}H'^{m.^). 
Then, multiplying (I2.4p i and (12.40 9 by ^1 and ^2 respectively, adding the resulting equations, and 
utilizing (l2.4P /i. we find that vr can be expressed by ■?/', i.e., 

TT = 7r(^,X3) 

= [-A/i(ei¥. + m" + (A'p + Mi? + M'iimv + m - M^eiiei V]/(Aien 

= [Xi^r - (A'p + A/.|eP + M'ilW - M^Giei V]/(A|er)- (2.40) 

Thus (I2.4p i can be rewritten as 

-if" + aip = uj (2.41) 

with boundary conditions 

ip{~oo) = ip{+oo) = 0, (2.42) 

where a = (AV+A/x|eP+M2e?)/(A/.) > ando; = ei[A/i^"' - (AV + A/i|eP + M'eiW]/iXfi\^\') 
G H°°{M.). By the ODE theory on a bounded interval and the domain expansion technique, we 
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can obtain a unique solution ip := ^pi^^x^i) E H'^{R) to (l^iTT) . ([232]). In view of ([230]), thus vr 
can be uniquely determined by the known functions ip and ip. Employing arguments similar to 
those in the construction of ip, we can obtain a unique solution 9 := 9{^, X3) G H°°{R) to (12 ■4p ^ 



with 9{—oo) = 9{+oo) = 0, where the solution 9 depends on the known functions n and ip. 

Finally, by a simple computation, we can check that {ip,9,ijj,n) with A constructed as above 
indeed is a solution to ([23]), <^M- ^ 

Remark 2.1. For each 0:3, it is easy to see that the solution {(p{C,, ■), 6'(^, ■), ip{^, ■), 7r(,^, ■), A(,^)) 
constructed in Theorem 12.21 possesses the following properties: 

(1) A(^), ip{^, ■) and 7r(^, ■) are even on ^1 or ^2, when another variable is fixed; 

(2) (p{^, ■) is odd on ^1, but even on ^2, when another variable is fixed; 

(3) 9{^, ■) is even on ^1, but odd on ^2, when another variable is fixed. 

Remark 2.2. We mention that the system (12. 4p . (12. 5 p for the vertical case enjoys rotational 
structure. Hence, we can also use the rotation method as in [16] to construct a solution of (12. 4p . 
(12.51) . which is simper than the above construction process. 

The next lemma provides an estimate for the if'^(R)-norm of the solution (yj, 6',?/', vr) with 
1^ varying, which will be useful in the next subsection. To emphasize the dependence on ^, we 
write these solutions as (</?, 6',?/', 7r)(^) := (</?, 6',-?/', 7r)(,^, X3). 

Lemma 2.1. Let k > 0, and D C A^ 6e a nonempty bounded set satisfying the closure D C A^. 
Then for any ^ G D there are positive constants A/^, B^, Ck and D^, which may depend on g, M, 
IX, p and D, such that 

mO\\HHm<^k, (2.43) 

MOWH^m^Bk, (2.44) 

||^(e)||H^(R)<Cfc, (2.45) 

lk(OllH^(R)<^fc, (2.46) 

where {ip,9,ip,n){^) and A(^) be constructed as in Theorem \2.i\ Moreover, 

II^(OIIW)>0. (2.47) 

Proof. We still give the proof for the horizonal case only. Throughout this proof, we denote 
by c a generic positive constant which may vary from line to line, and may depend on g, M, p, 
p and D. 

We begin with the estimate (I2.43p . We first rewrite (12. 6p as 



r\i) = [(A^P + 2A/i|er + (M6)')^"(0 + A'pV (0 

- \^\\X'p + A/i|er + (M^i)' - ^p')V^(0]/A/i, 



(2.48) 



which yields 



||^"(0||L^(M)=(A/i)-^ / [A(Ap + 2/x|er+(M6)V(0 + AVV^'(0 
- \^\'{X'p + Xp\e + {Mil? - gp')^{OW)dxs 



(2.49) 
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and 

IIV^'"(OIIl^w = -(Am)-^ / [A(Ap + 2/i|eP + (M6)^)V^"(0 + AV^'(0 

JR (2.50) 

- mX'p + A/i|eP + (M6)^ - 9p'M0]ri0dxs. 
Noting that since tp E A, the inequahty fl2.47p holds, and there is a constant c, such that 



HH 



< c. (2.51) 



On the other hand, in view of Proposition 12. 8[ we have 

A(0 > c > for any ^ e ©. (2.52) 

Thus, using Cauchy-Schwarz's inequahty, (12.5 ip and (I2.52p . we deduce from (12.490 and (I2.50p 
that 

whence, by (1^^ . 

\m{R) < c. (2.53) 



Consequently, differentiating (I2.48P with respect to xs and using (I2.53p . we find, by induction on 
k, that (12:431) holds for any A; > 0. 

Now we turn to the estimate of (12.440 . From (12.410 . we have 

Jr 
which, together with Cauchy-Schwarz's inequality and (12.430 . yields 

\MhHr) < c- 

We further deduce from (12.410 that 

ll<^ll?/2(K) < c. (2.54) 

Thus, we can employ (I2.54p and (I2.4ip to deduce that (I2.44p holds for any A; > 0. Using this 
fact, the estimates (12.43P and ( I2.44p . the expression of vr in (I2.40p implies (12.460 . Finally, similar 
to the deduction of ( 12.440 . we can show that ( I2.45P holds. This completes the proof. D 

2.4- Exponential growth rate 

In this subsection we will construct a linear real-valued solution to the linearized problem 
( 11.70 along with ( 11.60 which grows in-time in the Sobolev space of order k. 

Theorem 2.3. Let 

A= sup A(0, (2.55) 

?eAgn(L-iz)2 

then there exist a positive constant A* G (A/2, A), and a real-valued solution (^, u, N, g) to the 
linearized problem ([L7| ) along with U.6\) . such that 

(1) For every fc G N, 

l|(^,u,N,g)(0)||^. <oo (2.56) 
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(2) For every t > 0, {g, u, N, q)(t) G H'' , and 

e^**||(^,N,g)(0)||H. = ||(f5,N,g)(t)||^. = e*'^l(^,N,g)(0)||H., (2.57) 

e''*'\\ui{0)\\H, = \\uiit)\\H,=e'''*\\ui{0)\\H^, ^ = 1,2,3. (2.58) 

(3) In addition, 



divu(O) =divN(0) = 0, 



and 



11^3(0)1^2 >0. 
Remark 2.3. In view of the third conclusion in Proposition 12.21 and (I2.3ip . we see that 



(2.59) 
(2.60) 



0<A<supA(O<^Al|pVp|| 



Proof. Let 

v(^, xs) = -i(p{^, X3)ei - ie{^, 0:3)62 + V(^, 2:3)63, 

where {(p,6,ip) with an associated growth rate A(^) is constructed in Theorem 12.21 for any given 
^ e A^. Recalhng the definition of supremum, and using Proposition 12. 2[ we find that there exist 
^^ and ^^ := — ^^ such that 

A* := A(^^) = A(^2) ^ (A/2, A) and C G A^ for i = 1 and 2. 



In view of Remark 12.11 



f git,^) = -e^*VELi^3(r,a:3)e''"«'", 
u(t,x) = A*e^**ELi v(r,a:3)e-'«'", 
g(t,x) = A*e^**ELi^(r,a:3)e^'''«'", 

[ N(t, x) = e^** ELi M ■ Vx(v(r , a;3)e^'''«'") 



(2.61) 



gives a real- value solution to (ll.5p - (ll.7p satisfying (I2.57p - (l2.59p . Lemma \2A\ immediately implies 
that (^, u, N,g) constructed in (I2.6ip also satisfies (I2.56p . Finally, thanks to (I2.47p . we get 

||iZ3(0)||i. =2A* U\i\x,)dx, [ cos2(x'-e^)dx'>0. 

JR J{27rLT)2 



Hence (I2.60p holds. This completes the proof of Theorem 12. 3[ 
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Remark 2.4. When M = Me3, recalling the expression of M, one gets from a simple compu- 
tation that 



iNsimh = e'^*'\ 



N^ml^ = ^M'e'^'' I {d.,^{e,xs)fdxs / cos2(x' ■ e)d^' ■ 

J{2nLTf 



On the other hand, thanks to (12.470 and tp ^ H'^{R), we have dx^^ipi^^yXs) ^ 0. Consequently, 

||A^3(t)l|L2 = e''**||Ar3(0)|U2>0 foranyt>0. (2.62) 



3. Sharp growth rate of the hnear solutions 

In this section, we will show that A defined by fl2.55p is the sharp growth rate of arbitrary 
solutions to the linearized problem (11. 7p . Since the density varies, the spectrum of the linear 
operator is difficult to analyze, and it is hard to obtain the largest growth rate of the solution 



operator in some Sobolev space in the usual way. Instead, motivated by [10[, we use energy 
estimates to show that e^* is the sharp growth rate. First, we have the following energy identify. 



Lemma 3.1. Let (^, u, N,g) solve the problem i\1.6\) . ([i.Tp with M = Mj (i = l,2j. Then, 



— / (p|utp + M^l^iul^ - gp'ul) dx + 2/i / I V^iulMx = 0. (3.1) 

dt Jq Jn 

Proof. We differentiate (II. 7p in time, multiply the resulting equation by u^ and then integrate 
by part over fl. Using the first and third equations in (11.71) . we obtain (13. ip . D 

The next lemma allows us to estimate the energy in terms of A. 
Lemma 3.2. Let a real-value function u G H^ with divu = 0, then 

/ (^P>3p - M^l^iUp) dx < AM p|u|Mx + Ai2 iVupdx for i = 1,2. (3.2) 

Jn Jn Jn 

Proof. Since C'=^{n)nH^{n) is dense in H\n), it suffices to show (E^D for u G C°°{n)r\H^{n). 
Let / be the horizontal periodic Fourier transform of /, i.e., 

/(e, xs) = [ /(x', x,)e-'-'<d^' for ^ G (L-'Zy, 

J(27rLT)2 

and 

""ll^, X3) = -i^i^, Xs), U2(^, 3:3) = -iOi^, Xs), U3(^, X3) = ^(^, X3). 

Then ^icp + ^,2^ + ip' = because of divu = 0. Moreover, 

^lif iiO iiii) 
Vu = {diUj) = I 6V5 66' «6V^ 
-iif' -iO' Ip' 



By the Fubini and Parseval theorems, we have (see [10|, Section 3.3]) 

/^P>3rdx=^-^ Y^ gip'l{p'>0}+p'l{p'<0})\u3\'^dx3 

Jq ^ ge(L-iz)2 •^'^ 

^^ Yl / ^(p'l{p-'>o}+p'l{p-'<o})|V^|'dx3 (3.3) 



47r2L2 

CG(L-1Z)2 



21^ J2 ^p'l^rdxs, 



An 

5e(L-iz) 



where l{p'>o} and l{p'<o} denote the characteristic functions. Obviously, 

lim / gp'\ijj\^dx3 = / gp'\ij{0,X3)\'^dx3. 
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(3.4) 



In addition, we can control the following terms involved with ip. 



47r2L2 



47r2L2 



< 



47r2L2 



^ E /4(ien^r+ieiv^+6^ndx3 
E /leiriM^+i^r+i^ndxs 

3 

E E / l-^jl'd^s = 2 /" l^iupdx, 

dxs + E / 



ii. 



(ieriv^r+iv^'r)dx3 



47r2L2 



Ce(L-iz)2 



47r2L2 



?G{L-iZ)2 i= 



lim inf 



2L/,/|2 I U///|2^ 



(len^r + i^' 



2L/,/|2 I U/,//|2\ 



(ieii^r+i^' 



< 



47r2L2 



|t|2 •^ ' /L^ / |t|2 

'^' 5e(L-iz)2\{o}"^'^ '^' 

E /dv^'l' + l^f + l^'Hd^s = 2 / |93Updx, 



-da;3 



1 



47r2L2 



?e(L-iz) 



lim sup 



2L/,|2 I L/,/|2^ 



pimw + \i^ 



dX3 + E / 



2L/,|2 I L//|2^ 



Pderiv^p + iv^ 



< 



47r2L2 



'^' ^e(L-iz)2\{o}"^'* '^' 

E / PdV^I' + IV'I' + l^nd^3 = 2 / plupdx, 



-dx3 



and 



Ce(L-iz) 
1 



47r2L2 



I5K0 Jm \i\ 



+ E 

CG{L-1Z)2\{0} 






//|2^ 



-da;3 



47r2L2 



hmsup / m^i±mi±wf,^^ 



l?Ko 



+ E 

?G(L-1Z)2\{0} 



leP 



-da;3 



< 



47r2L2 



47r2L2 



E E / l^f da;3 = 2 /" iVupdx. 



5e(L-iZ) 



Ce(L-iZ)2 l<i,j<3' 



Now, let 



^(V^,0= / ^p'|^rda;3-M 



l?P 



/MeiMlV^P + l^ da;3, forM = Mei; 



A (l^? + ^) dx3, for M = Me3. 



(3.5) 



(3.6) 



(3.7) 



(3.^ 



(3.9) 
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By splitting Z{ip,$,) = ZiDl^ip),^) + Z(3{ip),$,), we may reduce the boundednes of Z{ijj,^) to 
the boundedness of Z {^{{ip) , ^) , since the imaginary part can be dealt with in the same manner. 
Hence, without loss of generality, assume that ip is real-valued. 

For ^ e AS n {L-^Zf y^ 0, we make use of (^M) and fl233D with s = A(0 to deduce that 



^(^>o<^ [pi\^m'+m')dx, + ^ I {Am^T + m'^ + rn dx. 



<- 



A 

w 



lel^ 



(3.10) 



. Jr Jr 

On the other hand, by Proposition 12. 3[ we infer that for any ^ ^ A^ \ {0}, 

A f pm'm' + m')dxs + fx [im'm' + m'^+rr)dxs 



^(V',0<o<^ 



A 



(3.11) 



Consequently, putting (I3.3l) - (l3.1ip together, we can conclude that (13.21) holds. In fact, taking 
i = 1 for example, we can estimate that 



2 [ {gp'\us\' - M'\d,u\') dx < -^ Y. I 9P'\M''dx, 



47r2L2 

1 

47r2L2 



lim inf 



lim sup 



^P 



'l„;.|2 






m'm' + w\') 



dX3+ J] Z(^,0 

?e(L-iz)2\{0} 



< 



1 



4.^LM"S!o''M5 



+ E 

?e(L-iz)2\{o} 



A 



A / pi\^m' + m')dx,+p im'm'+m'^+r\')dx, 

Jr Jr 

A / pm'm' + \^f)dxs + p /(4ienv^f + iieiv+v^'?)dx3 



<2A2 / p|u|Mx + 2A/i / iVupdx, 
which gives (13.21) for i = 1. 



n 



Now we are in a position to prove our main result of this section. 



Proposition 3.1. Let (^, u, N, g) solve (LI) with M = Mj (i = l,2j. Then we have the 
following estimates for any t > 0. 



f?(t)||i<Ce2^*(||f^o|li + ||uo||^. + ||VNo| 



L^j 



X = L'^ or H\ 



(3.12) 



|u(t)|| 



Ut 



tl'';ilL2 



+ r l|Vu(s)||i.ds < Ce'^'WlgoWh + l|uo||^2 + ||VNo||i.), (3-13) 

^0 



l|N(t)|U2 < Ce^"(||f?o||L2 + ||uo||h2 + llNoll^i; 
where the constant C depends on p, p, A and Mj. 



(3.14) 
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Proof. Integrating f l3.ip in time from to t, we get by Lemma [3.21 that 

\\y^ntit)\\l, + 2fi f \\VdMs)\\hds<Io + A'\\VP^mL2+A^i\\Vu{t)\\l, (3.15) 

Jo 

with 

Jo = ||A/pUi|i=o|li2 + \\MdiUo\\l2 -5f||M3(0)||^2, 

where the first term on the right-hand side of f l3.15p can bounded in terms of the initial data 
(^o,Uo,No). In fact, from (11.71) 9 and (11.21) we get 

/ p|u^|^(r)dx = / [(V X N) X M + /lAu - gge^] ■ u^(r)dx 
Ja Jn 

<C{Mr)\\h + l|Au(r)||i. + ||VN(r)||i.) + ^||u.(r)||i„ 
which implies 

\\Mr)\\h < C{Mt)\\1. + ||Au(r)||i. + ||VN(r)||i.) for any r > 0. 
Let r — > 0, we see that 

||ui|t=o|li2 < C(||f?o|li2 + ||Auo||i2 + llVNolli^). 
Integrating in time and using Cauchy-Schwarz's inequality, we find that 

All Vu(t) 11^2 =A||Vuo||i2 +2A f f Vdsu{s) : Vu(s)dxds 

Jo Jn 

<A||Vuo||i2+ / \\\/dsu{s)\\lAs + A^ [ ||Vu(s)||i2ds. 
Jo Jo 

On the other hand, 

Aa,||^/;^u(t)||i2 =2A /pu-a,u(t)dx< ||v^9,u(t)||i2+Alv^u(t)||i2. (3.17) 

Jn 

Hence, putting (I3.15p -f l3.17p together, we obtain the differential inequality 

d,\\^n{t)\\l. + fiWVumh <h + 2A (^l^umh + f^ 1^ ||Vu(.)||i2d.^ 
with Ji = 2/i||Vuo||L2 + Jo/A. An application of Gronwall's inequality then shows that 

\\V-Mt)\\h + f^ 1^ l|Vu(.)||i2ds <e^^iv^uo||i2 + ^ (e^^* - l) ^^^^^^ 

<Ce2^*(bo|li2 + ||uo||l,2 + ||VNo||i2) 
for any t > 0. Thus, making use of (13.151) . (I3.16P and (13.181) . we deduce that 

^IIV^u,(t)||i2+/.||Vu(t)||i2</i + A||v/;^u(t)||i2+2A/.y" ||Vu(5)||i2ds 

<Ce2^*(||f?o|li2 + ||uo||^2 + ||VNo||i2), 



(3.16) 
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which yields 03.131] . 

To prove (13.121) . we use f l3.13p and (ll.7P i to arrive at 



|^(t)||x <||^o||x 



/ ||^s(s)||xds 
Jo 

<||^o||x + ||p'||l°° / ||^3(s)||xds 



(3.19) 



<C^e^*(||^o|U + lluolk^ + IIVN0IU2) 
where X = L^ or H^. Similar to the derivation of (13.191) . one obtains 



|N(t)|U2<||No|U2+ / ||N,(s)|U2ds 
Jo 

<||No|U2 + ||M||loo / ||Vu(s)||L2ds 
Jo 



<C^e^*(||^o||L2 + ||uo||i/2 + ||No||i/i). 
Hence (I3.13P and (I3.14p hold. This completes the proof of the lemma. D 

4. Nonlinear energy estimates of the perturbed problem 

In this section, we derive some nonlinear energy estimates for the perturbed Cauchy problem 
(ll.4p - (ll.6p . and the integrand form of Gronwall's inequality of the high-order energy estimates. 
To this end, let {g,u,q) be a classical solution of the perturbed problem (ll.4p -( iL6l) with p : = 
^ + p > in (0, T) X f2 for some T > 0. Moreover, we assume that the classical solution {g, u, q) 
possesses proper regularity, which makes the procedure of formal deduction sense. 

In what follows, we denote by C{6o) a generic positive constant depending on 60. The notation 
a ^ b means that a < Ch for a universal constant C > 0, which may depend on some physical 
parameters in ([HD. 8 := £{t) := ^(^,u,N) := (||(p, N)||2^3 + Ml^Y^^ and Sq = ^(^o,Uo,No). 
a = (ai, 0^2, as) denotes a multi-index of order |a| = ai + 0^2 + 03. 

4.I. Estimate of the perturbation density 

We first note that by the classical Sobolev embedding results (see [l|. Chapter 5]), we have 

iklU4< IkllillVMllf, < ||m||hi, (4.1) 

\ML^<\MlM\h<\\u\\m, (4.2) 

Mw<M\Wh\\t?+^. (4.3) 

With the help of the above estimates, we can bound the perturbation density g. In fact, applying 
(9" to (ll.4P i. multiplying the resulting identity by d'^g in L^{^), we get 

Ijt E /(5"^)'dx = - J2 /5"(A3)5>dx+ Yl /5"(u-v^)a>dx 

0<|al<3"^^ 0<|a|<3"^^ 0<la|<3'^'^ (4.4) 



23 



where /i can be bounded as follows, using f l4.3p . Holder's and Young's inequalities. 

0<|q|<3 
^II^I|l2|I%IIl2 + ||V^||L2(||u3||i2 + ||VM3||l2) 

(4.5) 



I V'^|U2(||m3||l2 + IIVM3IU2 + \\V\3\\l2) 



+ \\V'g\\L2[\\u3\\L2 + IIVM3IU2 + \\V'us\\l2 + IIVXIU2) 

i i 11 

SIqWh^IWsWl^ + ||V^||i/2||n3||22||^3|lH2 + ||V^^||//1||U3||22||U3|||^3 

+ \\^^q\\l4u3\\12\\u3\\1, < CiSo)\\us\\l2 + So£^. 

On the other hand, keeping in mind that divu = 0, we use fl4.ip . f l4.2p and Holder's inequality 
to control I2 as follows. 

h <||Vu||i4||V^|U4||V^||L2 + ||V'u|U4||V^|U4||V'^|U2 + \\Vu\\Loo\\V^g\\l2 
+ ||V=^u||i4||Vf?|U4||V3f?|U2 + \\VMl^\\V^q\\l4^^0\\l^ 

+ ||Vu||Loo||V^^||i2+ Yl /u-a"V^a"^dx 

0<|a|<3"^^ (4.6) 

<||Vu||hi||V^||j:^i||V^||l2 + ||Vu||h2||V^||hi||V'^||l2 + ||Vu||H2||V'^||i2 



|Vu||h3||V^||^2-J V / |9>|Mivudx < £l 
2 n^i :^„Jn 



0<\a\<3 

Hence, one gets from (I4.4p - (l4.6p that 

1 d 



^^^.At)rH^<C{So)\\u,ml^+SoS'{t)+S%t). (4.7) 

In addition, we can deduce from 01.41) 1 that 

\\Qt\\L^<\\u3\\L2+£\ (4.8) 

||Vft|U2 < IIM3IU2 + II VM3IU2 +£^< C(5o)||%||l2 + So£ + S\ (4.9) 

\\Qu\\L^<\\^it\\Al + S)+S^ + S\ (4.10) 

4.2. Estimate of the perturbation magnetic field 

We continue to bound the perturbation magnetic field N. Applying 9" to ([I3D3! multiplying 
the resulting identity by (?"N in L^(r2), we obtain 

14 y /i9"N|2dx 

2dt ^ Ja ' 

= J2 /" 5"(M ■ Vu) • 9"Ndx + J2 [ ^"(^ ^ (" ^ ^)) ■ ^"Ndx 
+ ^ /" 9"(M ■ Vu) ■ (9"Ndx := J1 + J2 + J3, 

101=3"^^ 

where we have used the fact V x (u x M) = M • Vu. 
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Noting that V x (u x N) = N ■ Vu — u • VN, similar to the derivation of fl4.5p and f l4.6p . we 
can control Ji and J2 as follows. 



Jl 


< 


>: 


|5"VU|,.2| 


9"N| 


L2 














0<|a|<2 




















<\ 

r-^ 1 


Vu ^2 


|N| 


|l2+ I 


V^u| 


L^\ 


VN|i2 


+ |V=^U| ^2 


1 V^N|| 


L2 




<| 


1 1 

1 9 1 19 

U|l2|u||,, 


:|N|U. 


1 

> -1- 111 1 3 
' + ]"■ Il2 


2 

In 1 3 1 

|U|^3| 


|VN|i2 + | 


u\l,\u 


3 

1 4 
■1 i/4 




s 


0{So)\n\l. 


+ '^0^" 


5 













IV^NI 



L2 



and 



^2 <||Vu||i4||N||i4||N||i2 + (||Vu||i4||VN||i4 + ||V^u||L4||N||i4)||VN||i2 

+ (||Vu||ioo||V2N|U2 + ||V2u|U4||VN||i4 + ||V3u||i2||N|Uo.)||V2N||i2 

StvtII - I IIV72,,1| .1IV72t 



+ (||Vu|Uoo||V^N|U2 + ||V^u||i4||V^N|U4 

+ ||V^u|L2||VN||loo + ||V^u|L2||N||z,oo)||V^N|U2 



V / (u ■ a" VN) ■ a^Ndx 

^1 1/,"'^ 



<||Vu||hi||N||hi||N|L2 + ||Vu||H2||N||H2||VN||i2 + ||Vu||j:^2||N||^2||V^N 



0<|a|<3' 

il/l ||1-N ||i|/l||i"N 11^2 -h II VU||j|/2||i'N||j:/2|| ViN||L2 -t- II VU||j:/2||i\l||j:^2|| V i"N||L2 

+ ||Vu||H3||N||H3||V^N|L2-i V /" |9"NPdivudx < f I 
II WM II \\H II 11^ 2 ^^ /o 

Hence, one gets 

^|l|N(t)||^3 < C{6o)\\u{t)\\l,+6oS\t)+S%t) + J3. (4.11) 

In addition, we infer from fll.4P /L that 

l|N<|U2 <|| Vu|U2 +S^< C(5o)||u|U2 + 5oS + £'', (4.12) 

l|VN,|U2 < \\VM\l-+£' < Ci5o)\\u\\L2+6o£ + £^ (4.13) 

\\^tt\\L2 < ||uj^i(l + £)+£^ + £\ (4.14) 

^.5". Estimate of the perturbation velocity 

First we restrict the density p = f? + p to possess a positive lower bound. By virtue of (14. 2p . 
there is a constant ^q G (0, 1), such that 

inf^3eK{p(x3)} , 

II^oIIl- < ^ for any ||^o||h2 < 6^. (4.15) 

Consequently, in view of (ll.4p i and divu = 0, we find that for any (t,x) G [0,T) x fi, 

inf^3 6ffi{p(x3)} ^+ M ^ r+ ^ ^ r ^+ ^. ^ 3sup gR{p(x3)} 

< mf {po(t,x)} < p(t,x) < sup{po(t,x)} < '-- , 

I xen xen 2 

where po = ^0 + P- From now on, we always assume that £ < 5q < 5'q < 1. We remark that the 
upper and lower boundedness of p will be repeatedly used below. 
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We now estimate the time-derivative of the perturbation velocity. Multiplying f ll.4p 9 by u in 
L^(ri), utilizing fll.4p i and integrating by parts, we obtain 

--7- / (^ + p)|u|Mx + /i / |Vupdx= / [(Vx N) X (N + M)-^^e3] ■ udx 

Yq ■ udx, 



where Yq := (V x N) x (N + M) — gge^. Applying dt to fll.4p 9. multiplying the resulting identity 
by Ut in L'^{fi), we make use of fll.4P i. and integrate by parts to deduce 

7:-n / (^ + p)|ui|Mx + /i / |Vui|Mx 
^ o-t Jn Jn 

(V X Nj) X (N + M) + (V X N) X Ni - ^u ■ Vu 



- {g + p)ut ■ Vu - ggtes - gtUt 



Uidx 



Yi ■ Uidx. 



Applying d^ to fll.4p 9 again, multiplying the resulting identity by uf, employing fll.4P i and 



(V X N) X (N + M) = (N + M) ■ VN - iv|(N + M)\^, 



we have 



9QttG3 - QttUt 



ol7 / (^ + P)|uitpdx + /i / |Vuti|Mx 
2 dt J^ Jq 

= /" { - N,, ■ Vu,, ■ N - N, ■ Vu,, ■ N, - (N + M) ■ Vu„ ■ N„ 

- QtUtt - Qttu ■ Vu - gtUt ■ Vu - ^,u ■ Vu, - {g + p)u, ■ Vu, - {g + p)u„ ■ Vu • u„| 

2 

:= - V / 5;(N + M) ■ Vu„ ■ 9,2-^Ndx + / Y2 • u„dx. 
.•_n ^n ^Q 



dx 



Adding the above three equalities, we get 



2 2 

^1E /(^ + p)|5HMx + pX: / |V9iu|Mx 

j=0 '^'^ i=0 "^^ 

J2 [ ^i-dludx-J2 [ diCN + M)- Vu„ ■ 9,2-^ Ndx. 



(4.16) 



Obviously, 



|Yo||L2<||f?||L2 + ||N||^,||N 



L2||-'^ll//2 



INI 



< 



C(5o)||(f5,N)|U. + 5o^^ 



(4.17) 



Recalling that ^ < 5 < 1, we use dH]), (H^D, dlSD, dUD, fH:T2D and IKm to arrive at 

I|Yi||l2 <||^,||l2 + ||ft||L2||u||Loc||Vu||L°o + ||^ + p||Loo||u,||i2||Vu||Loo 

+ (1 + ||N|Uoc)||VN,||i2 + ||VN||i4||N,|U4 + ||ft|U4||u,|L4 



^ll^tlU^ 



u 



fl-3 



l^t||L2 + ||ft|ki||u,||Hi + ||Vu||^2||u,||i2 



(4.18) 



+ ||VN,|U2 + ||N||h2||N,||hi 
<C(5o)||u|U2+<5o^ + ^||u,||hi. 
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To bound Y2, we argue in a way similar to that used for (14.181) with the additional estimates 
fHTTU]) and d^D to infer that 



_L O 

IIY2IIL2 <||^tt||L2(l +£ + \\vit\\m) + ||ft||Hi(||ufi|||2||uii|||^i +£\\\it\\m) 
+ ||uf||/fi||Vui||//i +£:||uii||L2 
<(||ui||^i(l +£) + £^){l + £^ + IIu^IIhO + £{\\ntt\\l4viu\\\, + ^||ui||^2). 

Finally, using fl4T2|) -f H?T4|) . we obtain 

V / 9;(N + M) • Vuit ■ 9j2-^Ndx 

<[(l + £)||N«|U2 + ||Ni||^,]||Vuii|U2<[(l + ^)||ui||^i+£2]||u,,||^i. 

Putting fl4.16l) - fl4.20l) together and using Cauchy-Schwarz's inequality, we conclude 

2 2 



(4.19) 



(4.20) 



1 d 



9;^E(^ + /^)ll«^)lli^ + Ell^^*" 



2dt 



., < 



i=0 



j=0 



(4.21) 



+ (C(5o)||u(t)|U2 + 5o^(t) + £\t) + f (t)||u,(t)||H0l|ut(t)||L2 
+ [||u,(t)||^i(l + ^(t))+^2(t)](l+£2(t) + ||u,(t)||^2)||u,,(t)|U2 

+ £:(t)(||u,,(t)||i||u,,(t)||i+£(t)||u,(t)bOI|u**(t)|U2 

+ [{l + £{mMm\m+S\t)f. 
We proceed to estimate higher derivatives of the perturbation velocity. Rewriting (11.41) 9 as 
- /iAu +{q + p)ui + Vg = (V X N) X (N + M) - gge^ - {g + p)u ■ Vu, (4.22) 

and multiplying (14.221) by u^ in ^^(fi), we have 
/i d 



|Vu(t)|Mx+ / (f? + p)|ui|Mx= / [Yo - (f? + p)u ■ Vu] ■ utdx. 
2 "t Jf7 Jn Jn 

Differentiating (I4.22p with respect to t and multiplying the resulting equations by uu in ^^(fi), 
we obtain 

^- I \Vut{t)\M^+ I (g + p)\uu\^d^ = / [Yi - (^ + p)u ■ Vui] ■ utidx. 
Adding the above two equahties, we get 

f |E / iv9Mt)rdx+^ A,+p-)i9;+^updx 

i=0 '^'^ 1=0 "^^ 

= [Yo-{g + p)u ■ Vu] ■ Uidx + / [Yi - (^ + p)u • Vut] ■ Uiidx. 
Jn Jq 

Clearly, 



(4.23) 



[g + p)u ■ Vu]|i2 + ||(^ + p)u ■ Vui||i2 < £' + ^11 Vui]U2. 



(4.24) 
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Consequently, in view of fl4.18p . f l4.19p and fl4.24p . using Cauchy-Schwarz's inequality, we get 
from K^ that 

<C(5o)||(f.,u,N)(t)||i. + 5o^2(t)+^2(t)||u,(t)|||.. 



(-\^ 



Next, we continue to derive more estimates of higher-order derivatives of the perturbation 
velocity. Multiplying ( 11.4^ 9 by u^, integrating and using (14.3p . one gets 






2 , . c2 (4-26) 



while applying di to (11.40 9. multiplying the resulting equations by SjUj in L^(fi), and using (14. 3p . 
we have 



l|5«u,||i.<||y^+^9,u,||i.<C(5o)||(^,u,N)||i. + <5o^2 + (l + Ol|u*iiz.. 



2 

(4.27) 



<C(5o)||(f.,u,N)||i.+5o^l 
If we take did^ to (11.40 9 and multiply the resulting equations with diW-t in L^(r2), we see that 

In particular, summing up the above three estimates, we conclude 

\\M]i^<£'- (4.28) 

Finally, adding dOT]) to (I12SD, and utilizing (i2SD-(112HD, we arrive at 



\ 1=0 1=0 

1 



(4.29) 



+ ^||a^(Vu,u„Vu,)(t)||i.<C(5o)||(f.,u,N)(t)||i. + 5o^2(t), 



i=0 

provided 5o is sufficiently small. 

4.4- Energy estimates 

Now, we sum up the previous estimates (14.71) . (14. lip and (I4.29p . and use Cauchy-Schwarz's 
inequality to find that 



d 
It 



[q. N)(t)||^3 + Y^ ||V9;u(t)||i. + Y^ \\^ITpdlu{t)\\l. 

(4.30) 

+ ^||ai(Vu,u„Vu,)(t)||i. <C(<5o)||(^,u,N)(t)l|i. + J3 + <5o^2(t), 



1=0 

provided 5q is sufficiently small. 
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To deal with the term J3, we shall make use of the momentum equations fll.4p 9. Let a := 
(«!, 02, cts) be given and satisfy |a| = 3. Without loss of generality, assume a\ 7^ 0. By employing 
a partial integration, one sees 



9"(M- Vu) •(9°Ndx= / M ■ V(9^N ■ (9^udx, (4.31) 

'Q Jo. 

where 7 := {a\ — 1,02,03) and /3 := (ai + 1,02,03) satisfying I7I = 2 and |/3| = 4. Letting 
d^ := didj, and applying d'^ to f ll.4p 9. keeping in mind that 

(V X N) X M = M ■ VN - V(N ■ M), 

we find that 

{g + p)d'^ut + Wd'^q + gd'^ge^ — fiAd'^u 
= -utd^'ig + p) - di{g + p)djMt - dj{g + p)9iUj - d'<[{g + p)u ■ Vu] (4.32) 

+ a^[(V X N) X N] - V(a^N • M) + M ■ V(9^N. 

Multiplying f l02|) by 9^u with |/3| = 4 in L2(fi), we get 
-^ /"(^ + p)|9"u|2dx + /i /" |V9"uPdx 

= - [ di{g + p)d^ut ■ 9"udx + f utd\g + p) ■ 9^udx 

+ / di{g + p)djUt ■ d'^ud:x.+ / 9j(^ + p)9jUt ■ 9^udx 
in in 

+ / d^[{g + p)u ■ Vu] ■ 9'^udx - / 9^[(V x N) x N] ■ 9^udx 
Jn Jn 

9 

-g I d^gd^usd^ + [ {g + p)u- Vd"u ■ 9"udx - / M ■ V9^N • 9^udx = ^ /„ 
Jn Jn Jn ^^^ 

Recalling that £{t) < 60 <t^ 1, it is easy to verify that 



E'" 



n=l 



<\\dy\\L4d°us\\L^ + {\\d"vi\\L2 + £')\\uty2 + £' 



Therefore, 



1 d 
2di 



[g + pWu{t)\M^ + p / |V9"u(t)|Mx 



<C{6o)\\uml,+6oS'{t) + h, 

where we have used (14.281) in the last inequality. 
Adding fH:33|l to fH:30D and using fOTj) . we have 

1 2 



(4.33) 



d 
dt 



[^,N)(t)||^3+ 5^ \\V^+pd"u{t)\\l, + Y,\\'^dMt)\\l, + J2\\V^+pdMt)\\l. 



I«l=3 



i=0 



i=0 



+ 5]||aKVu,u„Vu,)(t)||i. + ||VMt)lli^<C^(^o)||(^,u,N)(t)||i. + 5o^2(t). 



i=0 
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In particular, integrating the above inequality over (0,t), we get immediately 



||(^,N)(t)|||3 + ||u(t)||i. + ||V^u 



2 
L2 



U/ 



U+ 



2 
L2 



5^ l|a:(vu,u„ vu,)(.)||i. + iivM^^"' 



L2 



j=0 



ds 



< ||(f?o,No)||^3 + ||uo||^3 + \\dtUofH2 + ll^^uollia 



(4.34) 



+ 



C(5o)||(^,u,N)(s)||i. + 5o^2^s) 



ds. 



Applying the classical regularity theory on the Stokes equations to f ll.4p 9. one sees that 

|V'ui||i2 + II Vgt|U2 <\\dt [(V X N) X (N + M) - gge^ - {g + p)u ■ Vu - (^ + p)ut] 1^2 

<(1 + ^')||ft|U2 + {S+ ||f?i|U2)||uJ^i + (1 + £)||Ni||^i + ||u,,|U2 (4.35) 
<£ + \\ut\\m + ||utJL2, 



and 

II V2u||h2 + II Vg||H2 <||(V X N) X (N + M) - gge^^ - {g + p)u ■Vu-{g + p)ut\\H2 

<£ + \\ut\\H2. 

On the other hand, letting t — )■ in (14.281) . we infer that 



\\dtno\\l.<S, 



I q2 ||2 < c2 



Similarly, it's easy to infer that 

'% L10IIL2 ^ <^0- 

Consequently, combining (I4.34p with (I4.35p -( l4.38p . we conclude 

£2(t) + ||ui(t)||^2 + ||Vgi||i2 + ||Vg||^2 + WuttWh 



(4.36) 

(4.37) 
(4.38) 



+ 



< £^ 







1 



\h 



ds 



J2\\dl{Vu,u,,Vu,){s)\\l, + \\VMs) 

i=0 

/ [C(5o)||(^) u, N)(s)||^2 + 6o£'^{s)~\ ds, provided Sq is sufficiently small. 
Jo 



In addition, from (14.281) we get the following estimate, which will be used in Section [51 



IL2 



+ IIGI 



H2 



Hl|2 < C2 
\\m ^ ^ ' 



where F = -u- V^, G := (V x N) x N - (^ + p)u • Vu - gUt, and H := V x (u x N). 

Now, we mention that the local existence of classical solutions and global existence of classical 
small solutions to the 3D nonhomogeneous incompressible Navier-Stokes equations (i.e. M = 
in (11.11) ) have been established by many authors, see 29|, |32(| for example. However, the 
mathematical analysis becomes much more difficult for the MHD equations (ll.ip due to the 
loss of magnetic diffusivity (dissipation). To our best knowledge, there are only two existence 
results on the 3D compressible MHD problem ([Hi]) i- ([HI]) 3 without gravitational field. The 
first one is on the local existence of classical solutions by Kawashima |19, Theorem 6.14], while 
the second one is on the local existence of strong solutions to (ll.ip i- Ql.ip q by Li, et al. |2J, 
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Theorem 1.1]. By a proper modification, one can adopt basic ideas of tlie fixed point theory 
as in J2J, Theorem 1.1] to show that there exists a unique local classical solution (^, u, N) G 
C{[0,T],H^{n) X H^{n) xH^{n)) to the perturbed Cauchy problem ([Ll])-([L6D for any given 
initial data satisfying £^(^0,^0, Nq) < cxd. Moreover, the classical solution satisfies the estimates 
above if £{t) < 60 on (0,T). The proof of the local existence on a horizontal periodic domain is 



standard (please refer to (9|, |ll| for well-posedness results of other models on horizontal periodic 



domains), in fact, it is easier than that for the corresponding compressible case, hence we omit 
it here. Finally, summing up the above estimates and the existence statement, we arrive at the 
following conclusion: 

Proposition 4.1. For any given initial data {qq, Uq, Nq) G H^ x H^ x H^ satisfying f l4.15p and 
divuo = divNo = 0, there exist aT > and a classical unique solution {g, u, N) G C([0, T], H^ x 
H'^ X H^) to the perturbed problem f ll.4p - fll.6p . Moreover, there is a constant 60 G (0,1], such 
that if S{t) < 6q on [0,T], then the classical solution satisfies 

£\t) + \\ntit)\\j,, + ||(Vg, nt)t\\h + l|Vg||^2 



+ / \j2\mVu,u,,Vu,){s)\\l. + \\VMs)\\l^ 



ds 



i=0 



(4.39) 



<Ci6o)£^+ [ [C((5o)||(f?,u,N)(s) 
Jo 



lL2 



+ A*£'(s) ds 



and 



\\F\\l, + \\G\\l, + M\l,<C{5o)£', 
where A* is provided by Theorem \2.3[ 



(4.40) 



5. Proof of Theorem 11.11 

Now we are in a position to prove Theorem II. II by adopting the ideas in [8|, [l5|, l3l|- First, in 
view of Theorem 12.31 we can construct a linear solution (^', u', N') G H°° to (ll.7p with the initial 
data (^0)Uo,No) G H°° and divuo = divNg = 0; moreover, the linear solution satisfies 



.A*t 



Uo 





("■] 




= e^*' 


( '° \ 


VNo/ 


c 


\w) 


c 




VNo/ 



A* G (A/2, A), 



(5.1) 



< e 



A*t 



u\\\l2 and ||mo3||l2 = 1, 



(5.2) 



where the sharp linear growth rate A > and the linear norm || ■ ||£ are defined by ( I2.55P and 




1^1112 + 



u 



m 



INIP 



respectively, and M30 and u\ stand for the third component of uq and u', respectively. 



Denote Cq 

No), and Ci := 



[^o,Uo,No)||£: 



[^o,Uo,No)^| 



s := ^(^-o,Uo,No), (^S,ug,N^ 



^(^0,Uo, 



sical solution [g^ ,\i^ , N*^ 



f)o, Uo, No)||l2. By virtue of Proposition 14.11 there exists a unique local clas- 
G C([0,T],if3 X H'^ X H^) to flLi|) . emanating from the initial data 



[gi ui N^) with 



U,<^i) 



le = CqS. Let us now choose 60 G (0, 1) as small as in Proposition 
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14.11 and let C{So) > be the constant appearing in Proposition 14. II for the fixed choice of S^. Let 
L = min{(5o, ^o}; and 6 G (0, l) satisfy 



T^= lln^>0, i.e. 5e^*^* = 2e, 
A* 6 



0) 



(5.3) 



where the value of Eq will be fixed later to be sufficiently small and independent of 6. We then 
define 

u' I it) < 6o 



T* = sup{te (0, T" 



and 



T** = sup^tG (0,T" 



N* 



u' |(t) 



< 26Cie 



A*t 



L2 



where T™^^ denotes the maximal time of existence. Obviously, T*T** > 0, and furthermore, 



u^ I [T* 



6o and 






25Cie 



AT* 



(5.4) 



L2 



if T* < oo and T** < cx). Then for all t < min{r^ T*, T**}, we deduce from the estimate IK39\\ . 
and the definitions of T* and T** that 



Q 

u' \{t) 



<C{6o)6' 



+ C{6o 



No 



<A* 



<A* 



u"^ I (sj 

N'5 



A* 



u"^ I [S] 

2 




2 

L2 



u*^ I (s: 



ds 



ds 



2 r2^2A*t 



ds + CoC{5o)5' + AC{5o)Ct5'e 



ds + Ca^'e 



2„2A*t 



(5.5) 



for some constant C2 := C2(5o) > independent of S. Then applying Gronwall's inequality to 
the above estimate, one obtains 



u' |(t) 



io 



(5.6) 



for any t < min{T^T*,T**}. 
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Let (^^, u'^, N'^) = {g^, u^, N^) - S{g\ u\ N'). Then {g^ u^ N;^) := S{g\ u\ N^) is also a linear 
solution to (11.71) with the initial data (^q,Uq,Nq) G H°°. Moreover, (f?*^, u*^, N*^) satisfies the 
following non-homogenous equations: 



r ^d I ^/„,d 



gt + p'u1 = F\ 

puf + Vq'^ + gg'^e3 - /iAu"^ - (V x N)<^ x M = G^ 

Nf - V X (u'i X M) = H^ 

divu^ = 0, divN'i = 0, 



where 



F' 



-n' . Vg' 

(V X N"^) X N-^ - {g^ + p)u^ ■ Vu"^ - ^u^^ 
V X fu-^ X N-^) 



with initial conditions 

(/(O), u'i(O), N'i(O)) = 0, divu^ = and divN^j = 0. 
From the estimates (14.401) and (15. 6p . it follows that 



F' 
G' 

a' 



< 2C2C{So)6'e''^ * for any t < min{^^ T*,T**}. 



Now, with the help of Proposition 13. II and Duhamel's principle, we follow a procedure similar 
to that in [7|, Section 4] to infer that 



g'^ 

u'i |(t) 



Jo 



L2 



G^ 



ds 



(5.7) 



2C,C{6o)6' 2A** 



Jo 



2A*-A 



where we have used the fact that 2 A* — A > 0. 
Now, we claim that 

T'' = min{T^T*,T**}, 

provided that small Eq is taken to be 

_ . / 5o Ci 1 



(5. 



(5.9) 



Indeed, if T* = mm{T\T*,T**}, then by ([5SD and (ESD we have 



Q 
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which contradicts with fl5.4p . On the other hand, if T** = mm{T^,T*,T**}, in view of 05. ip . 
3]), ( 1^ and the fact that A/2 < A* < A, we have 




L2 



+ C-ad e 



A*T** 



(Ci + 2CsEo) 



<26Cie 

which also contradicts with (15 ■4p . Therefore, f l5.8p holds. 
Finally, we again use fl5.9p . fl5.7p and f l5.2p to deduce that 



\ui{T')\\,.>\\uUT')h^-\\ut{T 



L2 



6\\ul{T')\\,. - \\ul{T')\\,. 



(5.10) 



where ujjg and U3(T'') denote the third component of ujj and u'^(T'^), respectively. This completes 
the proof of Theorem II. 1[ 

Acknowledgements. The research of Fei Jiang was supported by the NSFC (Grant No. 
11101044 and 11271051), the research of Song Jiang by NSFC (Grant No. 11229101) and the 
National Basic Research Program under the Grant 2011CB309705. 



References 



[1 

[2: 
[3 
K 

[6 

[r 

[8 
[9 

[10: 

[11 

[12: 



R. A. Adams, J.F.Fournier John, Sobolev Space, 2005, Academic Press: New York. 

H. Cabannes, Theoretical Magnetofluiddynamics, Academic Press, New York, 1970. 

T.G. CowHng, Magnetohydrodynamics, Interscience Publishers, New York, 1957. 

R. Duan, F. Jiang, S. Jiang, Rayleigh- Taylor instability for compressible rotating flows, submitted, 

2011. 

D. Erban, The equations of motion of a perfect fluid with free boundary are not well posed, Comm. 

PDF 12 (1987), 1175-1201. 

Y. Guo, C. Hallstrom, D. Spirn, Dynamics near unstable, interfacial fluids, Comm. Math. Phys. 

270 (2007), 635-689. 

Y. Guo, H. Hwang, On the dynamical Rayleigh- Taylor instability. Arch. Rational Mech. Anal. 167 

(2003), 235-253. 

Y. Guo, W. Strauss, Instability of periodic BGK equilibria, Comm. Pure Appl. Math. 48 (1995), 

861-894. 

Y. Guo, I. Tice, Almost exponential decay of periodic viscous surface waves without surface tension. 

Arch. Rational Mech. Anal. 207 (2013), 459-531. 

Y. Guo, I. Tice, Linear Rayleigh- Taylor instability for viscous, compressible fluids, SIAM J. Math. 

Anal. 42 (2011), 1688-1720. 

Y. Guo, I. Tice, Local well-posedness of the viscous surface wave problem without surface tension. 

Anal PDF 2012, to appear. 

Y. Guo, I. Tice, Compressible, inviscid Rayleigh- Taylor instability, Indiana Univ. Math. J. 60 

(2011), 677-712. 



34 



[13] R. Hide, Waves in a heavy, viscous, incompressible, electrically conducting fluid of variable density, 

in the presence of a magnetic field, Proc. Roy. Soc. (London) A 233 (1955), 376-396. 
[14] H. Hwang, Variational approach to nonlinear gravity-driven instability in a MHD setting. Quart. 

Appl. Math. 66 (2008), 303-324. 
[15] J. Jang, I. Tice, Instability theory of the Navier-Stokes-Poisson equations, 

arXiv:1105.5128v2[math.AP] 13 Jun 2011 (2011). 
[16] F. Jiang, S. Jiang, G. Ni, Nonlinear instability for nonhomogeneous incompressible viscous fluids. 

Science China Math. 56 (2013), 665-686. 
[17] F. Jiang, S. Jiang, W. Wang, On the Rayleigh- Taylor instability for the incompressible viscous 

magnetohydrodynamic equations, arXiv:1204.6402v2 [math.GM] 1 May 2012. 
[18] F. Jiang, S. Jiang, W. Wang, On the Rayleigh- Taylor instability for two uniform viscous incom- 
pressible flows, submitted, 2012. 
[19] S. Kawashima, Systems of a hyperbolic-parabolic composite type, with applications to the equations 

of magnetohydrodynamics, Ph. D. Thesis, Kyoto University, 1983. 
[20] M. Kruskal, M. Schwarzschild, Some instabilities of a completely ionized plasma, Proc. Roy. Soc. 

(London) A 233 (1954), 348-360. 
[21] A.G. Kulikovskiy, G.A. Lyubimov, Magnetohydrodynamics, Addison-Wesley, Reading, MA, 1965. 
[22] J. Larmor, How could a rotating body such as the sun become a magnet. Rep. Brit. Assoc Adv. 

Sci, 159, 1919. 
[23] L.D. Landau, E.M. Lifshitz, L.P. Pitaevskii, Electrodynamics of Continuous Media, Vol.8, 1984 

(Translated from the Russian). 
[24] X. Li, N. Su, D. Wang, Local strong solution to the compressible magnetohydrodynmic flow with 

large data, J. Hyper. Diff. Eqns. 8 (2011), 4315-4369. 
[25] F. Lin, P. Zhang, Global small solutions to MHD type system (I): 3-D case, Comm. Pure. Appl. 

Math. 154 (2013) (to appear). 
[26] L. Nirenberg, On elliptic partial differential equations, Estratto dagli Ann. del. Scuola Norm. Super. 

Pisa Ser. HI, XIII, Fasc. II, 1959. 
[27] J. Priiss, G. Simonett, On the Rayleigh- Taylor instability for the two-phase Navier-Stokes equations, 

Indiana Univ. Math. J. 59 (2010), 1853-1871. 
[28] L. Rayleigh, Analytic solutions of the Rayleigh equations for linear density profiles, Proc. London. 

Math. Soc. 14 (1883), 170-177. 
[29] L. Tong, H. Yuan, Classical solutions to Navier-Stokes equations for nonhomogeneous incompress- 
ible fluids with non-negative densities, J. Math. Anal. Appl. 362 (2010), 476-504. 
[30] Y. Wang, Critical magnetic number in the MHD Rayleigh- Taylor instability, J. Math. Phys. 53 

(2012), 073701. 
[31] Y. Wang, I. Tice, The viscous surface-internal wave problem: nonlinear Rayleigh- Taylor instability, 

Comm. PDE 37 (2012), 1967-2028. 
[32] P. Zhang, Global smooth solutions to the 2D nonhomogeneous Navier-Stokes equations, Inter. 

Math. Research Notices, Article ID rnn098, doi:10.1093/imrn/rnn098, 2008. 



35 



